Bound states in the one-dimensional two-particle Hubbard model with an impurity 
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We investigate bound states in the one-dimensional two-particle Bose-Hubbard model with an 
attractive {V > 0) impurity potential. This is a one-dimensional, discrete analogy of the hydrogen 
negative ion H~ problem. There are several different types of bound states in this system, each of 
which appears in a specific region. For given V, there exists a (positive) critical value Ud of U, 
below which the ground state is a bound state. Interestingly, close to the critical value {U < Uci), 
the ground state can be described by the Chandrasekhar-type variational wave function, which was 
initially proposed for H~. For U > Uci, the ground state is no longer a bound state. However, 
there exists a second (larger) critical value Uc2 of U, above which a molecule-type bound state is 
established and stabilized by the repulsion. We have also tried to solve for the eigenstates of the 
model using the Bethe ansatz. The model possesses a global Z2-symmetry (parity) which allows 
classification of all eigenstates into even and odd ones. It is found that all states with odd-parity 
have the Bethe form, but none of the states in the even-parity sector. This allows us to identify 
analytically two odd-parity bound states, which appear in the parameter regions —2V < U < —V 
and —V < U < 0, respectively. Remarkably, the latter one can be embedded in the continuum 
spectrum with appropriate parameters. Moreover, in part of these regions, there exists an even- 
parity bound state accompanying the corresponding odd-parity bound state with almost the same 
energy. 

PACS numbers: 03.65.Ge, 03.75.-b, 71.10.Fd 



I. INTRODUCTION 

The hydrogen atom serves as a paradigm in quantum 
physics. The exact solution 1] of the two-body problem 
shows that one proton can bind one electron. However, 
a very interesting and highly nontrivial problem is the 
following: How many electrons can be added, if at all, to 
the hydrogen atom to form negatively charged ions? The 
answer, thanks to the endeavor of many people during 
more than half a century 04^, is: at most one, yield- 
ing the fascinating object H~ (for some very readable 
reviews, see 0, Q)- The difficulty of this problem lies 
in the fact that the Coulomb repulsion between the elec- 
trons is of the same order of magnitude as the Coulomb 
attraction of the proton, which makes it a delicate job 
to balance the kinetic energy and the repulsion against 
the attraction to allow a bound state with two electrons 
below the continuum threshold. Actually, the repulsion 
between the electrons is so crucial that H~ has only a 
single bound state i.e. the ground state, in which the 
two electrons are strongly correlated 

The problem can be generalized to any other atom: 
How many electrons can bind to a given nucleus of charge 
+Z7 As a fundamental problem in both atomic and 
mathematical physics, it has attracted a lot of interest 
and many profound results have been established Q . 

With the advent of the era of cold atoms, it becomes 
natural to pose similar questions for cold atoms in opti- 
cal lattices. It is well-known that in a one-dimensional 
tight binding model, a defect potential — < creates 
a localized mode around the defect (10|. In the absence 



of atom-atom interaction, an arbitrary number of bosons 
can be trapped in this localized mode. However, as soon 
as an on-site atom-atom repulsion J7 > is turned on, at 
most a finite number of bosons can be trapped simulta- 
neously, since the interaction energy scales with N"^ while 
the single particle energy scales with N, if N bosons re- 
side in the defect mode. The problem is then the follow- 
ing: What is the maximal value of N for given values 
of V and U7 A different but related problem concerns 
the existence of bound states for fixed N and varying V 
and U. In particular one may ask whether the ground 
state is a bound state. Intuitively speaking, there must 
be a critical value of U beyond which the ground state 
becomes delocalized by the repulsion. 

These problems motivated us to study the bound states 
in a Bose-Hubbard model with an attractive impurity po- 
tential, for which some results were reported in Ref. [ill : 
in the limit of large U, V a similar model was considered 
recently in Ref. 12,. Since numerical calculation is indis- 
pensable and since the dimension of the Hilbert space 
scales as with M being the number of sites in a fi- 
nite lattice, we shall only consider the second problem 
above and the case N = 2. Therefore, we have a one- 
dimensional, discrete analogy of the H~ problem. The 
attractive defect potential plays the role of the central 
Coulomb attraction while the atom-atom interaction cor- 
responds to the Coulomb repulsion between the electrons. 
However, the possible experimental realizations with cold 
atoms allow more freedom in the choice of V and U, 
which are freely adjustable. Especially, both of them can 
be either positive or negative. 

It turns out that this simple model embodies rich 
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physics. For example, unlike the electron-electron repul- 
sion in H~, here the repulsion between the atoms plays 
a dual role. On the one hand, it may destroy a bound 
ground state, on the other hand, it allows formation of 
some other bound state. Specifically, for a given attrac- 
tive potential {V > 0), there exists a critical value Ud 
of U, above which the bound ground state is destroyed 
by the atom-atom repulsion. But for U > Uc2 > Uci, a 
molecule-type bound state appears. This is a repulsion- 
aided bound state — the two atoms are bound into a 
molecule by the repulsion and are trapped as a whole 
by the impurity potential. A noteworthy feature of this 
bound state is that it can be identified as a hound state at 
threshold (T3l - [l7| . which is normalizable (thus of a finite 
size) even at the continuum threshold. 

A most unexpected object realized in this model is a 
"bound state in the continuum" (BIG) [11]. It is a nor- 
malizable state yet its energy is located within the con- 
tinuum spectrum of the model. For N — 2, the model 
seems simple enough to try a diagonalization via the 
Bethe ansatz , although the continuum version of the 
model is known to be nonintegrable (20l |. It turns out 
that in our case a Z2-symmetry (parity) divides the total 
Hilbert space into two sectors, one of which (odd un- 
der parity) allows diagonalization with the Bethe ansatz, 
while the other exhibits diffractive behavior. This en- 
ables us to construct the two odd-parity bound states in 
closed form, one of which can be a BIG. 

This paper is organized as follows. In Sec. |lll the 
model is defined and some of its basic properties are dis- 
cussed. In Sec. mil the numerical algorithm for identi- 
fying bound states in the system is introduced and the 
typical bound states found are displayed. Sec. IIVI and 
IVl treat the repulsive delocalization of the ground state 
and the repulsive localization of the molecule-type bound 
state, respectively. In particular, the values of Uci,2 are 
determined both numerically and variationally, and the 
critical behavior of the bound states is investigated and 
compared. In Sec. IVI[ we discuss the Bethe ansatz solu- 
tion of the model and the two analytically solvable odd- 
parity bound states. The accompanying bound states to 
the two analytically obtained bound states are the sub- 
ject of Sec. IVIII We summarize our results in Sec. IVlIIl 



II. THE TWO-PARTICLE BOSE-HUBBARD 
MODEL WITH A DEFECT 



We consider two (spinless) interacting bosons (or 
equivalently, two spin- 1/2 fermions in the spin singlet 
space) in a one-dimensional lattice with a site defect. 
The Hamiltonian of the system reads. 
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FIG. 1. (Color online) The three continuum bands revealed 
by using the two quantities oi D = \xi\ + \x2\ and Pout- 
The three horizontal lines are clearly visible in each panel. 
The eigenstates and eigenvalues En are obtained by exact di- 
agonalization on a (21/ + l)-site {L = 20) lattice and with 
open boundary conditions. The values of the parameters are 
{V, U) = (1.5, 5.5), and R = L/2. 



Here cix (aj.) is the annihilation (creation) operator for a 
boson at site x. The Hamiltonian is just the conven- 
tional Bose-Hubbard model but with an impurity po- 
tential at site 0. The value of the impurity potential 
is ~V while the on-site interaction between two bosons 
is U. Note that the hopping strength is set to unity. 
The Hamiltonian is invariant under the refiection about 



the defect (da;,aj,) — >■ (a_ 



which defines the parity 
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of each state. Furthermore, the canonical transformation 
(a^;, aj.) — > (— )^(a2;, aj.) relates the eigenvalues and eigen- 
states of the Hamiltonian with parameters (V, U) to those 
of the Hamiltonian with parameters {—V,—U). There- 
fore, we shall consider in the following only the attractive 
impurity case > 0. The Hilbert space is spanned by 
Fock states which are defined as \xi,X2) = (>'ti'^x2\^) 
xi < X2, or \xi,X2) = (a^i)^ |0) if a^i = X2- 

The Hamiltonian ([T]) is written in second quantized 
form. For purposes below, it is also useful to work 
in the first quantization formalism. In this formalism, 
the wave function will be denoted as f{xi,X2) with xi^2 
being two independent variables. The Bose symmetry 
requires that f{xi,X2) — f{x2jXi). The correspond- 
ing wave function in second quantization is of the form 
V' = T,xi<x2 iPixi,X2)\xi,X2), with 

i^{xi,X2)^ [\/2-(\/2-l)4i,.J /(xi,X2). (2) 

The Hamiltonian is defined by its action on the wave 
function 

Hf{xi,X2) ^ - ^ [f{xi+j,X2) + f{xi,X2+j)] 

i=±i 

+ [-V{S.x,,o + S^.^o) + US^.^x,] f{xi,X2). (3) 
A state f{xi,X2) is even/odd under parity if it satisfies 

f{xi,X2) = ±f{-Xi, -X2)- 
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A. Three continuum bands 



We aim to identify the bound states, i.e., states in 
which the two particles are localized around the impurity. 
However, as a reference, it is necessary to identify the 
three continuum bands in this system first. 

The first band corresponds to the situation that the 
two bosons are neither captured by the impurity nor 
bound together by the interaction between them. For 
this band, the impurity and interaction cause phase shifts 
but do not contribute to the energy of the system, which 
is given by the kinetic energy of the two bosons, and thus 
this band covers the interval [—4, +4]. 

For the second band, one boson is captured by the 
impurity yet the other is not. The energy of the first 
boson is -VV"2 +4 ^ 

while that of the other boson 
lies between —2 and 2. This band covers therefore the 
interval [-VV^ +4-2, ~VV^ +4 + 2]. 

The third band corresponds to a dclocalized molecule 
state [2l|: The two boson are bound together by the 
interaction between them and the composite moves as 
a whole on the lattice p2j. Again, the impurity causes 
phase shifts or local modifications on the wave function 
yet does not contribute to the energy. This band covers 
[-VTPTTG, [/] if [/ < or [U, + 16] if C/ > [2l|. 

The presence of the three bands can be demonstrated 
numerically by using some quantities which can reveal 
the distinct nature of the extended states in them. The 
first quantity is the sum of the distance of the two par- 
ticles to the defect, Z? = jxil + |x2|. Suppose we choose 
a lattice of 2L + 1 sites with L sites on each side of 
the defect. For the first and third band, since the par- 
ticles move through the lattice either independently or 
bound together, (D) should be ~ L. For the second 
band, since always one particle is localized around the 
defect, (D) is expected to be ~ L/2. Another quantity 
is the probability of finding at least one particle out- 
side a ball with radius R and centered at the defect, 

PoutU) = E„.ax{|.,M..|}>.?, 1/(^1. ^2)P. If i? = i/2, 

it is easy to see that Pout should be ~ 0.5 for the sec- 
ond and third band, while for the first band it should 
be ^ 0.75. These predictions are readily verified numer- 
ically. In Fig. [U we see how the three bands, although 
overlapping in energy, are separated by using D and Pout ■ 
Moreover, their band edges coincide with the predicted 
values very well. 

Here some remarks are necessary. In ((Ij and ([3]), the 
Hamiltonian is defined on an infinite lattice, because 
bound states and extended states are only then well de- 
fined. However, in the numerical simulations, we can only 
deal with finite-sized lattices. We shall always choose a 
finite lattice with M = 2L + 1 sites, open boundary con- 
ditions and the defect located in the middle, as in Fig. [T] 



The Hamiltonian will be denoted as Hm and reads 

X——L X——L 

-Valdo. (4) 

Note that the defect site is set in the middle to preserve 
the reflection symmetry of the model. Thus, all the eigen- 
states have definite parity. In particular, it is easy to 
prove that the ground state [written as |g) or fg{xi^X2)\ 
is even and positive everywhere up to a global phase. 



III. BOUND STATE IDENTIFICATION 

The quantities D and Pout can not only be used to 
identify the three different bands, but also to discern the 
bound states from the extended states ^23(| . The latter are 
the overwhelming majority. The point is that there is an 
essential difference between bound states and extended 
states in their response to the change of the boundaries. 
An extended state would be very sensitive to the size 
of the lattice or the conditions (periodic or anti-periodic 
or in-between) imposed at the ends of the lattice: Its 
energy should change continuously as the boundary con- 
ditions change continuously and the typical amplitude of 
the wave function should decrease polynomially as the 
lattice size increases. On the other hand, the energy of 
a bound state should be independent of the lattice size 
or the boundary conditions, at least if the lattice size is 
much larger than the characteristic size of the wave func- 
tion, since its wave function decays rapidly away from the 
center. In the case of extended states, the value of {D) 
would increase linearly proportional to M (as argued in 
the previous section), and the value of Pout (with a fixed 
R, however large it is) would converge to unity, as M 
increases to infinity. But for a bound state, the value of 
{D) would converge to a finite value and the value of Pout 
would converge to a value which is always smaller than 
unity and decrease with growing R. 

Thus, to identify a bound state, we use the follow- 
ing algorithm: Begin with a lattice of Mi = 2Li + 1 
sites with the defect site in the middle, and open bound- 
ary conditions. Solve all the eigenstates and eigenener- 
gies of Hmi numerically by exact diagonalization 
and calculate the corresponding values of D and Pout for 
each eigenstate. Denote the n-th eigenstate as 
its eigenenergy as E^'^\ and the corresponding values 
of D and Pout as D^^'^ and Pn^^\ respectively. Here 
1 < n < dim{Mi), with dim{Mi) = Mi{2Mi + 1) 
being the dimension of the Hilbert space on the Mi- 
sited lattice. The same procedure is then repeated on 
a (sufficiently) larger lattice of M2 = 2L2 + 1 sites, with 
L2 > Li. A bound state would be signified by a pair 
of fixed values of D and Pout, at some fixed eigenen- 
ergy. More precisely, there should exist two indexes 
ni and 712 associated with each bound state, such that 
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FIG. 2. (Color online) A bound state (indicated by the arrow 
in each panel, where a blue * and a red o coincide) found 
by using the lattice expansion algorithm. In each panel, the 
blue * markers correspond to the smaller lattice with 51 sites, 
while the red o markers to the larger lattice with 71 sites. The 
values of the parameters are {V, U) = (1.5, 7), and R — 12. 

(i?if \P(^)) ^ (i?(f^\i?if ),Pr^)), respec- 

lively. 

As an illustration, we have carried out this algorithm 
for a specific pair of values of (V, U) on a smaller lattice 
(Li = 25) and a larger lattice {L2 = 35). The resulting 
values of D and Pout versus the eigenenergies are plot- 
ted in the upper and lower panels of Fig. [51 respectively. 
There we see that as the lattice is enlarged, the markers 
in the three continuum bands move upwards significantly, 
in accord with our expectation. However, in each panel, 
in the gap between the first and third band and at the 
same energy, we have a blue marker coinciding with a red 
marker, which indicates a bound state. Actually, it is a 
repulsion-aided bound state, which will be discussed in 
detail in Sec. IVl below. 

Finally, we would like to mention that the idea of 
distinguishing localized states from extended states by 
their different sensitivities to boundary conditions was 
employed before in the context of Anderson localization 
[25]. Besides that, it was also proposed as a criterion 
of differentiating superfluid states from insulator states 
[26^ . A superfluid state would response to the twist of 
the boundary conditions by generating a finite current, 
while an insulator state would remain inert under this 
distortion. 



A. A list of typical bound states 

With the effectiveness of the algorithm proven, we have 
explored the (V,U) parameter space extensively. Some 
general behavior can be observed and is illustrated in 
Fig. [3] by taking a cut through the parameter space at 
fixed V — 1.5 and varying U. We can see how the bound 
states change as U is varied from large negative to large 
positive values. Only the quantity D is displayed, but 
the same conclusions are obtained by using Pout- 

(1) In Fig. [3^, we see that for sufficiently large negative 



U, the ground state is a bound state and the only one. 
The fact that the ground state is localized is expected, 
since the interaction between the bosons strengthens the 
binding effect of the attractive potential. 

(2) In Fig. we see that as U is increased to —2, 
there are now two more bound states besides the ground 
state — the first and second excited state (indicated by 
the arrow). Further investigation shows that the first 
excited state is always a bound state and is of odd-parity 
whenever —2V < U < —V. Actually, this bound state 
is analytically solvable by the Bethe ansatz. It will be 
discussed in detail in Sec. IVII As for the other bound 
state, detailed investigation shows that it appears only 
in a subset of the —2V < U < —V region. In a certain 
sense, it is a companion of the analytically solvable one, 
always of even-parity and slightly higher in energy. It 
will be discussed in detail in Sec. IVIII 

(3) As U is further increased to —1, as in Fig. [5};, the 
two (excited) bound states disappear and are replaced by 
another bound state (indicated by the arrow). This is no 
ordinary bound state since its energy falls in the [—4, -1-4] 
continuum band. It is a "bound state in the continuum" 
(BIG). Further numerical and analytic investigation re- 
veals that it appears whenever —V<U < 0, and has 
odd-parity and the Bethe form as the analytically solv- 
able bound state in Fig. [Sb- They will be discussed to- 
gether in Sec. |VT1 

(4) As U is increased to 2, as in Fig. the ground 
state is still a bound state, and the only one. However, 
its energy, which has been increasing since Fig. is 
now very close to the lower edge of the second band. 
Here we should mention the simple fact that in the 
non-interacting case t/ = 0, the ground state is the 
only bound state. It has the simple form fg{xi,X2) oc 
M^i)M x2), with (X e--/« (e = l/\n[^{V + 
vV^ + 4)]) being the single particle localized mode in- 
duced by the defect. 

(5) As C/ > grows, a point is reached where the 
ground state is no longer bound. Actually, for J7 = 5, 
shown in Fig. [3^, there is no bound state at all. We 
may conclude from Fig. [SJi to Fig. [3^ the existence of 
a critical value of U, denoted as Ud, beyond which the 
ground state is turned into an extended one. 

(6) However, as U becomes even larger, as in Fig. [3f 
with U = 6.5, a new bound state (indicated by the arrow) 
appears. It is located in the gap between the first and 
third band. In contrast to the ground state, this bound 
state is established and stabilized by the strong repulsion 
between the bosons. Obviously, there must be a second 
critical value Uc2 of U, above which this type of bound 
state becomes possible. 

IV. STABILITY OF THE GROUND STATE 

As expected and demonstrated in Fig. [31 a sufficiently 
large U can delocalize the ground state \g). In this sec- 
tion, we are interested in the critical value Uci as a func- 
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FIG. 3. (Color online) Typical bound states found by the filtering algorithm. From (a) to (f), the value of U increases 
monotonically while the value of V is fixed to 1.5. The sizes of the smaller (blue *) and larger (red o) lattice are Mi = 51 and 
M2 = 71, respectively. In each panel, a pair of blue * and red o which coincide with each other corresponds to a bound state. 
In (a), only the ground state is a bound state. In (b), besides the ground state, the first and second excited states (indicated 
by the arrow) are bound. In (c), besides the ground state, there is another bound state (indicated by the arrow) which is 
embedded in the continuum band. In (d), again only the ground state is bound. In (e), there is no bound state at all. In (f), 
a bound state (indicated by the arrow) appears in the gap between the first and third band. 



tion of V and the critical behavior of the ground state. 
First, we note that for a positive value of t/, the second 
band has the lowest lower band edge among all the three 
bands (see Fig. [T] as an example). Moreover, this up- 
per limit for the ground state ener gy depe nds only on V 
and it has the value Eedge2 = —VV^ + 4 — 2. In terms 
of energy, it is then expected that as U increases, the 
ground state energy increases according to the Hellmann- 
Feynman theorem, up to the lower band edge of the sec- 
ond band and finally, at [/ = Ud, it coalesces with it. 
This can already be seen in Fig. [3] and is most clear in 
Fig. HJi, where the spectral graph of the model is shown. 

More information is obtained by studying the evolu- 
tion of the the ground state wave function fg{xi,X2) as 
U increases. In Figs.|3^-|3]:, we have plotted the profile of 
fg{xi,X2) at three increasing values of U. In Fig.|4^, we 
see that for large negative U, the attractive defect poten- 
tial is reinforced by the attractive interaction and the two 
bosons are sharply localized at the origin. In Fig. HId, for 
zero interaction U = 0, fg is simply a product of the wave 
functions of the two bosons. The critical behavior of the 
ground state is exhibited in Fig. where U is close to 
Uci- There, fg takes on a cross shape, which means that 
one boson is tightly bound by the defect, while the other 
is much more loosely bound. This indicates that the 
ground state wave function becomes delocalized continu- 
ously, since the wave functions in the second band have 
the cross shape as well, except that they are infinitely 
extended along the axes. This behavior is quite reason- 
able physically: When one boson is localized around the 
defect, the other boson feels only a screened potential, 
and thus is relatively loosely bound. 

Interestingly, a similar scenario occurs also for the hy- 
drogen negative ion H~ . Due to the difficulty mentioned 
in the introduction, initially it took a great effort to prove 
H~ as a bound system. However, Chandrasekhar had the 
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FIG. 4. (Color online) (a)-(c): Images of the ground state 
wave function fg{xi,X2) at three increasing values of U. 
Specifically, U = —2, 0, and 1 from left to right, (d): the 
bottom part of the spectral graph of the model as U increases 
towards the first critical value Ud- The value of V is set to 
0.5. The lattice has 201 sites. In (d), the lowest blue line 
is the ground state level and the dotted horizontal line in- 
dicates the bottom threshold of the second continuum band 



insight that H~ should have an "in-out" structure such 
that whenever one electron is close to the nucleus, the 
other is kept far away. He thusproposed a simple two- 
parameter trial wave function the essential part of 
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which reads 

X(ri,r2) oc e-"'^!-^'-^ +e-"'-^-'5'-\ (5) 

where presumably a S> /3. It turns out that the energy 
minimum is attained at a = 1.03925 and (3 = 0.28309, 
and the corresponding energy is sufficient to prove bind- 
ing for . 

Motivated by the observation of the critical behavior 
of the ground state and the similarity of the model with 
, we have tried the following variational wave function 
for the ground state: 

f{xi,X2) oc (^e~l^ilAi-|^2|/A2 ^g-|xi|/A2-|r.2|/Ai^ 

X (l -se-l^i-^^l/^^') . (6) 

Here Ai^2.3 are some characteristic lengths and < s < 1 
is an attenuation factor. We take the convention of 
Ai < A2. Here, besides the Chandrasekhar part (the 
first line), an additional factor (the second line) is in- 
troduced to suppress the wave function when the two 
bosons are close to each other. This is motivated by a 
fine structure of the wave function in Fig. H) That is, 
the value of the wave function at {x, —x) is visibly larger 
than its value at {x^x). It turns out that the ground 
state can be described by the variational form ^ quite 
well. In Fig. [SI the ground states obtained by exact diag- 
onalization are compared with those obtained within the 
variational scheme of We see that they agree even 
in details. Quantitatively, the overlap between the exact 
and variational ground state exceeds 0.99 in both cases. 
Note that we did not treat s as an independent parameter 
in the variational calculation. Instead, it is determined 
as a function of A3, i.e., s = \/ {1 - e-'^/^^)/{l + e'^/As). 
However, the agreement is very good even without vari- 
ation of s. 

The critical value Ud can be determined efficiently by 
exact-diagonalization. We choose a M-sited lattice and 
using the Lanczos algorithm we can solve for the ground 
state energy Eg{U,V;M) quickly. By the Hellmann- 
Feynman theorem, this quantity is a monotonically in- 
creasing function of U. Therefore, the critical value Ud 
is uniquely determined by the condition Eg — Eedge2- 
By using the Newton bisection method, it can be pinned 
down quickly. The Ud determined in this way depends 
of course on M . Actually, since a finite lattice can be un- 
derstood as an infinite lattice with infinitely large poten- 
tial beyond some radius, Eg{U, V; M) is a monotonically 
decreasing function of M by the Hellman-Feynman the- 
orem. This means the Ud determined on a finite lattice 
is lower than the true value. However, for M sufficiently 
large, the error is inappreciable. We have also tried to 
determine the transition line variationally: If it is possi- 
ble to find for a given U a set of Ai_2,3 and s such that 
the energy of © is lower than Eedge2 , then U \s a, lower 
bound oiUd- Both the exact diagonalization result (solid 
line) and the variational one (dashed line) are shown in 
Fig. ini There we see that the two approaches agree with 
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FIG. 5. (Color online) (a) and (b): images of the ground 
state wave functions obtained by exact diagonalization on a 
finite lattice with 301 sites. The parameters are (a) [V, U) = 
(0.5,1) and (b) {V,U) = (1.5,2.5). (c) and (d): images of 
the variational wave function minimizing the energy, with the 
same parameters as in (a) and (b), respectively. The inner 
product of the exact wave function and the variational one is 
0.991 and 0.998, for (a)-(c) pair and (b)-(d) pair, respectively. 
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FIG. 6. (Color online) The bound-extended transition line 
of the ground state determined by various means. Blue solid 
line: exact-diagonalization on a 301-sited lattice; cyan dashed 
line: modified Chandrasekhar variational wave function (see 
Eq. ([SI, s 7^ 0); red dotted line: non- modified Chandrasekhar 
variational wave function (s = 0); green dash-dotted line: 
Uci = F as in App. \^ 



each other very well in the full region of V investigated. 
This again demonstrates that the wave function near the 
critical point can be well described by the variational 
form (|6]). 

The modified (or improved) Chandrasekhar wave func- 
tion ^ is accurate but not very tractable analytically. 
Therefore, we have also tried some simpler variational 
schemes. The results are less accurate but capture well 
the qualitative behavior of Ud as a function of V . First, 
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in the double-occupied subspace is then 





(a) V = 0.1 


(b) y = 0.2 


(c) y = 1.8 

n 





Xi 





Xi 







FIG. 7. (Color online) Images of the squared wave function 
fmoii^iy ^2) of the molecule-type bound state with three dif- 
ferent values of V while the value of U is fixed to 7. The wave 
functions are obtained by exact diagonalization on a lattice 
of 201 sites. In (a) and (b), the wave function decays expo- 
nentially along the diagonal xi = X2- 



by a Hartree approximation as in App. \^ we can prove 
that Uci > V. Second, by the non-modified Chan- 
drasekhar (i.e., s = 0) approximation as in App. [Bl we 
can show that Ud > 1.5V as V —?' and Uci > V + 1 
as V — )■ -1-00. These two resuhs are also displayed in 
Fig. [6] We see that they are not so accurate as ([6]) , but 
still capture the behavior of Ud to leading order in V 
as y — >■ +00. The fact that the non-modified Chan- 
drasekhar wave function is less accurate than the mod- 
ified one indicates that the suppression factor is quite 
essential near the critical point. 



REPULSION-AIDED BOUND STATE 



In the preceeding section, we discussed the repulsion 
driven delocalization of the ground state. However, as 
revealed in Fig. [3f, the repulsion can also facilitate the 
formation of a localized state. For a given V, there exists 
a critical U (denoted as Uc2), beyond which a bound 
state appears in the gap between the first and third band. 
The mechanism for its formation relies on the repulsion 
between the two bosons which can bind them together 
to form a molecule (possible in a tight- binding model), 
which moves on the lattice as a whole with some effective 
hopping strength, and can thus be trapped by the defect 
potential. The picture is most clear in the large U limit. 
In this limit, the double-occupied states \x)m = \x,x) 
are degenerate and higher in energy than the non-double- 
occupied states by a gap of order U. Due to transitions 
through intermediate states like \x,x + 1), two adjacent 
double-occupied states \x)m and \x + l)m are coupled by 
an effective hopping strength 2/U, and each of them is 
up-shifted in energy by A/U. The effective Hamiltonian 



+00 



^eff ^ ^ {\x)rn7n{x+l\ + \x + l),n,n{x\ 



U 



(7) 



which has the form of a standard tight-binding model. 
The dispersion relation reads E,n{k) — U + -^(1 + cosk), 
with a minimum at fc = tt. This molecule band is essen- 
tially the third band identified in Fig. [TJ Actually, its 
band range [U, U -\- 8/J7] is approximately the exact one 
[[/, vC/^-i-16] in the limit oiU ^ +00. Now turn on the 
defect potential. If the potential is weak enough, we may 
still confine the analysis to the double-occupied subspace. 
The defect is then equivalent to a term —2V\Q)mm{^\- As 
is well known [I^] , this leads to the formation of an even- 
parity bound state below the band minimum. Its energy 
is 



En 



U 



4 
U 



/4I/2 



2 



4 
U 

if^«|. 



(8) 



If [/ > 4, Emoi would be in the gap between the first 
and the third band (i.e. 4 < Emoi < U), for sufficiently 
small V. Thus, the bound state is expected to be stable. 
To verify the argument above, we show in Fig. [7^-b the 
squared wave functions of the bound state for two small 
values of V (and a fixed value of U). The wave function is 
significantly different from zero only on the line xi — xi 
and decays exponentially away from (0, 0), which are two 
salient features in accord with the molecule picture. 

In the analysis above, the defect potential is assumed 
to be weak enough to allow the restriction to the double- 
occupied subspace. \iV > 1, modifications are neces- 
sary. In this case, the ratio of the direct coupling (= v^) 
between the molecule state |0)„i and the disassociated 
states |0,±1) to the gap between them ([/ — F ~ U), is 
larger than the ratio of the effective coupling (~ 2/C/) 
between the molecule states |0)m and | ± !),„ to the gap 
between them (= 2V). Therefore, as a first approxima- 
tion, we assume that the statejO)m gets dressed by the 
state ^(|0,+1) + I - 1,0)) This is confirmed in 

Fig- Eti where V is increased to 1.8. We see that the 
wave function is no longer extended in the Xi = X2 direc- 
tion, but takes on a cross shape — the values of the wave 
function / at (0, ±1) and (±1,0) are much larger than at 
(1,1) or (—1,-1). Diagonalizing the Hamiltonian in the 
two dimensional subspace of {|0)m, -^(|0, +1) + | — 1, 0))}, 
we get the energy of the bound state 



1 



Ernol ^^{U-W) + v/(f/ - 1^)2/4 + 4 



U-2V 



U -V 



if J7 - V" > 4. 



(9) 
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FIG. 8. (Color online) The second critical value of (7 as a 
function of the defect potential strength V. The blue solid 
line is obtained numerically by using the lattice expansion 
algorithm, with a smaller lattice of 201 sites and a larger 
lattice of 301 sites. The green dashed (red dotted) line is the 
small- (large-) V approximation as in Eq. pop . Note that Uc2 
starts from 4. 




7.24 7.26 7.28 7.3 

U 



FIG. 9. (Color online) A portion of the spectrum (only even- 
parity levels are shown) of the model. The value of V is set to 
2. The lattice is of 101 sites with open boundary conditions. 
The highest one in the figure is the bound molecule one. The 
narrow anti-crossing is indicated by the arrow. Note the scales 
of the axes. 



Again, for sufficiently large U, Emoi lies in the gap be- 
tween the first and third band. 

By equating the value of Emoi to 4, we obtain the crit- 
ical value Uc2 for small and large V: 

Uc2^ < 4 (10) 

The value of Uc2 can also be determined numerically 
by using the lattice expansion algorithm. The result is 
shown in Fig. [5] (the solid line). We see that as V in- 
creases, the value of Uc2 interpolates smoothly between 
the small- and large- F behavior in Eq. (ITU)) , with the 
crossover occurring around ~ 1. These results confirm 
the pictures above. 




3 3.2 3.4 3.6 3.8 4 4.2 7.285 7.29 7.295 

U U 



FIG. 10. (Color online) The probability of finding both bosons 
inside a radius of J? = 10, (a) in the ground state and (b) 
in the molecule-type bound state. In both panels, different 
markers indicate different lattice sizes (blue *: 101 sites; red 
o: 201 sites; green x: 301 sites). The value of V is set to 2. 
The vertical dotted line in (a) indicates the critical value Ud ■ 
Note the scale of the horizontal axis in (b). 



A. A bound state at threshold 

We have shown that the repulsion can both destroy the 
ground state as a bound state and establish the molecule- 
type bound state. It is important to note that there is 
an essential difference between the two delocalization- 
localization transitions. The difference is, first of all, ex- 
hibited in the spectral graph. In Fig. [HJ the molecule 
bound state energy dX U ^ Uc2 is shown together with 
the few highest levels in the first band. In comparison 
with Fig.HJl, where the ground state level merges into the 
second band smoothly without appreciable anti-crossing, 
here the molecule level merges into the first band by a 
very narrow anti-crossing. The different behavior of the 
levels implies that the first transition is a continuous one, 
while the second transition is essentially a discontinuous 
one. We have seen in Sec. lIVI that during the first transi- 
tion the wave function of the ground state becomes more 
and more extended along the xi^2 axes as U increases 
towards Ud from below. Eventually, at U — Ud, the 
characteristic length A2 diverges and the ground state 
turns into an extended state residing in the second band. 
The whole process is continuous. However, for the sec- 
ond transition, as U decreases towards Uc2 from above, 
the wave function of the bound state is always of a finite 
(and small) size. Actually, the good agreement between 
the two-level prediction [Eq. ([TU)) ] and the exact value 
of Uc2 for large V indicates that as long as the bound 
state is there, it can be well approximated using the two- 
dimensional subspace above. It happens in a very narrow 
interval of U that the molecule bound state turns into an 
extended state, in which the two bosons are both delo- 
calized and move virtually independently. 

The two different pictures of the two transitions can be 
best demonstrated by using the quantity Pin = 1 — Pout, 
i.e. the probability of finding both bosons within a radius 
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R from the defect. In Fig. [TUb . Pin is plotted versus U 
for the ground state on three lattices with different sizes. 
We see that Pm is a continuous, monotonically decreasing 
function of U, and also a monotonically decreasing func- 
tion of the lattice size. Arguably, on an infinite lattice. 
Pin will drop continuously to at Ud- The behavior 
of Pin is thus reminiscent of the order parameter in a 
second order phase transition. In contrast, in Fig. llOb . 
where the value of Pin for the molecule-type bound state 
is plotted. Pin transits from almost vanishing to almost 
unity in a very narrow interval. Moreover, the transi- 
tion region gets smaller and smaller as the lattice size 
increases, which strongly indicates that the transition is 
discontinuous on an infinite lattice This behavior is anal- 
ogous to a first order phase transition. The discontinuity 
can be gleaned from Fig. [7] actually. There, for a given 
U > 4:, a.s V increases from 0+ (transversing the phase 
diagram of Fig. [5] from left to right), the wave function of 
the molecule-type bound state shrinks. This trend is in- 
compatible with the configuration of the wave functions 
in the first band, where the two atoms are both delo- 
calized and move virtually independently. Therefore, the 
transition between the two types of configurations cannot 
be continuous. 

It is an intriguing fact that a bound state becomes 
delocalized abruptly. Though so far we have failed to 
prove this rigorously, we note that similar scenarios were 
noted previously in various single-particle [Tsl - fTEj or few- 
particle [l^, [13 systems, and especially, in the H~ prob- 
lem (Trj . It actually occurs in a simple and familiar 
single-particle model — a three-dimensional finite depth 
square potential well ^is'] . Due to the repulsive core cre- 
ated by the centrifugal barrier, a bound state with zero 
energy or a bound state at threshold is possible 

in this model for angular momentum I > 1. The re- 
markable property of this type of state is that it sits on 
the threshold to the continuum yet has finite size. This 
means that if we adjust the parameter of the potential 
(e.g. the well depth) continuously across its critical value, 
we would observe a discontinuous behavior of P,;„ similar 
to Fig. [IU}d. 



VI. TWO ANALYTICALLY SOLVABLE BOUND 
STATES 

In Fig. [5};, we have a bound state inside the [—4, -1-4] 
continuum band. This is absolutely unexpected, since 
generally a bound state should lie outside of the con- 
tinuum spectrum. The reason is best demonstrated by 
Mott's argument for the existence of sharp mobility edges 
[2?! . In case of degeneracy between a localized state and 
an extended state, an infinitesimal perturbation would 
mix them and convert the localized state into an ex- 
tended one. It means that although such exotic objects 
called "bound state in the continuum (BIG)" do exist, 
as pointed out by von Neumann and Wigner as early 
as in 1929 [l^, they are expected to be unstable under 




FIG. 11. (Color online) The X1-X2 lattice decomposed into 
six regions by the a;i-axis, a;2-axis, and xi = X2 line. Note 
that adjacent regions share the boundary between them. Es- 
pecially, the origin belongs to all the six regions. 



perturbations. 

Nevertheless, it became recently possible to realize 
them experimentally in an electronic system psf . In 
our model, exploration of the {V, U) space shows that 
the BIG is quite robust. It is always present whenever 
—V<U < 0, though not necessarily inside the contin- 
uum. The robustness of the BIG indicates that it is pro- 
tected by a hidden, perhaps even strong symmetry. This 
motivates the attempt to solve the model analytically 
using the Bethe ansatz. It turns out that this approach 
allows us to understand both the BIG in Fig. 131:; and one 
of the bound states in Fig. |3]3. 

A. The Bethe ansatz approach 

We work with Eq. ([3]) and try to solve the eigenstates 
using the Bethe ansatz. An important observation is 
that the impurity potential and the interaction between 
the two particles are effective only on the three lines of 
2^1,2 = and xi — x^- Away from these lines, we have 
free particles. This observation suggests to try an ansatz 
for the eigenstates of the Bethe form which is character- 
ized by just two parameters, k\ and ki. Specifically, in 
region Ii (0 < cci < 2:2) in Fig. [Til the wave function is 
postulated to read 

/(a;i,X2)= 51 XI XI ^J'K.'Ti.a^) 

<To=0,l 0-1=04 (T2=0,l 

X exp [i(-)'^ifciXi+^o + i{c-y'^k2X2-a,^ , (H) 

where j((To, ci, (T2) = 4iTo + 2cri + <T2 + 1- The wave func- 
tion in regions I2 and I3 are defined similarly but with 
the A's replaced by i?'s and C's, respectively. The value 
of the wave function in the other three regions is then 
determined by the condition fi(x\^X2) — f{x2,xi). In 
each region, we have eight different plane waves. The 
reason is that the interaction between the two particles 
can exchange their momenta, and the impurity potential 
can reverse the momentum of a particle. 

With the wave function in the form above, the eigen- 
value equation Hf = Ef, with E = — 2cosfci — 2cosfc2, 



10 



is satisfied away from tfie three lines. Now we just need 
to get it satisfied also on the three lines. This requires, 
besides the eigenvalue equation above, first of all, that 
/ be single valued on the interface between two neigh- 
boring regions. Specifically, on the boundary between Ii 
and I2, we have 



and 



J+2 



B 



J+2) 



'V{B,+B, 



+2) 



(12a) 
1,2, (12b) 



and 



Aj + ^i+i ~ Bj + Bj+i, 
e^'^HA, - B,) + e-^'^HA,+, - B,+,) 



(13a) 



-V{B,+B,+,), J =5, 7. (13b) 



On the boundary between I2 and I3, we have 



Cj+C,+i=B,+B,+i, (14a) 
e^'=ni?.--Q) + e-''=^(S,+i-Q+i) 

^ ~V{Bj + Bj+,), J = 1,3, (14b) 



+ Cj+2 - Bj + Bj+2 , 



(15a) 



e*(S,-Q) + e-'^ni3.+2-Q+2) 



— iki I 



= -V{B,+B,+2), J = 5, 6. (15b) 

On the boundary between Ii and ![, and the boundary 
between I3 and I3, we have 

2i[(-)"sinfci + (-)"sinfc2](^j+4 - A,) 

= U{Aj + Aj+4), (16a) 

2i[(-)™sinfci + (-)"sinfc2](Cj+4 - Cj) 

= U{Cj+Cj+i), (16b) 

with 171,71 = 0,1 and j = 2m + n + 1 + (— )". Equa- 
tions (fT ^ - (fTB|) constitute a set of 24 homogeneous hnear 
equations of the 24 unknows ^'s, B's, and C's. For this 
set of equations to have nontrivial solutions, the 24 x 24 
coefficient matrix, which depends on ki, ^2, V, and U, 
should be singular. It is verified that the matrix is indeed 
always singular with a kernel of dimension 2. 

Here instead of dealing with the 24 x 24 matrix, some 
simplification is possible. Note that the system is reflec- 
tion invariant with respect to the impurity, which means 
if f{xi,X2) satisfies the eigenvalue equation Hf = Ef, 
then so does f{—xi, —2:2). Therefore, we can classify the 
eigenstates into even and odd ones. 

For the even case, f{xi,X2) = f{—xi,—X2), we need 
Bi ~ Bg-i and Ai = Cg-i. Now, it is ready to verify 
that if we have Bi = Bg^i, then Eqs. (fT^ - P^ lead to 
Ai — Cg-i automatically, and if Eq. (fTBb) is satisfied, 
then p6b ) is also satisfied. Solving the ^'s in terms of the 
B's from (fT2)) and (fT3)) . and substituting the expressions 
into P^ ). we get 



{U + 2i (sin fci - sin fc2 ) ) yBi - 
{U — 2z(sin/ci — sinfc2)) ^i?4 - 
{U + 2i(sin ki + sin fc2)) ^-82 - 
{U - 2i(sin ki + sin fc2)) ( B3 



V 



2i sin ki 
V 

2i sin fci 
V 

2i sin fci 
V 

2i sin fci 



(Bi + B3) 
{B2 + S4) 
{B2 + B4) 
(Bi + Bi 



+ {U - 2i(sinfci - sinfc2)) yB^ - 
+ {U + 2i (sin fci - sinfc2)) (^i - 
+ {U ~ 2i(sinfci -I- sinfc2)) - 
+ (J7 + 2i(sinfci +sinfc2)) { B2 



V 



2i sin fc2 
V 

2i sin fc2 
V 

2i sin fc2 
V 

2i sin fc2 



(S3 + B4) 
{B1+B2) 
{B3 + B4) 
{Bi + B2 



0, (17a) 
0,(17b) 
0, (17c) 
= 0.(17d) 



We thus get a 4 x 4 matrix, which should be singular 
to admit non-trivial solutions of Bi, B2, B^, B^. But the 
problem is that the matrix is always non-sigular as long 
as UV 7^ 0, since its determinant is —6AU^V^. Therefore, 
the even-parity eigenstates do not comply with the ansatz 

of inD. 

For the odd case, f(xi,X2) — — /(— xi, — 2:2), we need 



Bi — —Bg^i and Ai = — Cg_i. Similarly, it is ready to 
verify that if we have B^ — —Bg^i, then Eqs. ((T^ - 
lead to Ai = —Cg-i automatically, and if Eq. (fTOk ) is 
satisfied, then (fTBb 'l is also satisfied. Solving the A's in 
terms of the -B's from (|12p and (IT^ . and substituting the 
expressions into p^ ). we get 
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{U + 2i(sin ki ~ sin k^)) yB\ - 
{U ~ 2i{smki - sinfc2)) + 
{U + 2i(sin ki + sin fc2)) - 
{U - 2i(sin ki + sin fc2)) ( -B3 + 



V 



2i sin fci 
V 

2i sin fci 
V 

2i sin fci 
2i sin fci 



(Si + B3] 
(B2 + 
(B2 + Bi] 
{Bi + B3 



{U - 2i(sinfci - sinfca)) ^Bi - 
([/ + 2i(sinfci - sinfc2)) (^i + 
([/ - 2i (sin fci + sinfc2)) + 
- (L/ + 2i(sinfci +sinfc2)) { B2 - 



V 



This time, the 4x4 matrix is always singular and has a 
kernel of dimension 2. The value of 2 may be explained 
with the time reversal symmetry of the Hamiltonian: If / 
is an eigenstate, then so is /* with the same eigenvalue. 
This entails that for real fci. 2, if (i?i, S2, S3, B4) is a so- 
lution of (|T8)) . then {B\. B^, B2, B^) is also a solution of 

Here some remarks are necessary. Suppose we allow 
both Bose and Fermi symmetry. Then Eq. ^ itself is 
invariant under both the exchange of xi X2 and the re- 
flection ofa;i_2 — >■ — a;i,2. The two symmetries divide the 
Hilbert space into four sectors. For the anti-symmetric 
(fermionic) sectors, the interaction is ineffective and we 
have essentially free fermions in an impurity potential. 
The wave functions are in the Slater form and hence also 
in the Bcthc form. Therefore, only in one of the four sec- 
tors, i.e. the symmetric (bosonic) sector with even parity, 
the wave functions are diffractive and not in the Bethe 
form. Note that the presence of diffraction in a related 
model defined on a continuous line was shown almost 
fifty years ago in a classic paper by McGuire [1^. The 
interesting point here is that the reflection symmetry de- 
composes the bosonic Hilbert space into two subspaces, 
of which one shows no diffraction and can therefore be 
considered integrable. This is confirmed by an analysis 
of the algebra of scattering matrices and the associated 
Yang-Baxter relations, which is the subject of the next 
section. 



B. The integrable odd- parity subspace 

We shall now formulate the problem in terms of scat- 
tering theory, usually employed in the analysis of one- 
dimensional many-particle systems. This will explain 
the possibility to find eigenstates which have the Bethe 
form although the bosonic model is nonintegrable even in 
the two-particle sector. To this end we write the ansatz 
Eq. prj) in all regions as follows: 

/(a;i,a;2) = 

S Xfl.(2;i,^2X.,(A:i,fc2)e'('^^'=^^^+'^^'=^^^). (19) 

R,cri ,"'2 



2i sin fc2 
V 

2i sin fc2 
V 

2i sin fc2 
V 

2i sin fc2 



(S3 + B4) 
(S1+B2) 
{B3 + B4) 
[Bi + B2] 



= 0, (18a) 
0,(18b) 
= 0, (18c) 
= 0.(18d) 



I 

Here, Xr{^1t^2) are characteristic functions of the six 
regions Re{[i,j,k]}, 

[i, j, k] -ir^ Xi < Xj < Xk, i, J, fc G {0, 1, 2}, with xq — 0, 

while tJj e {1,-1} denote a pseudospin degree of free- 
dom accounting for the backscattering generated by the 
potential —V. Finally, the operator S symmetrizes the 
wavefunction with respect to xi and X2- The explicit 
use of S allows to reduce the number of amplitudes in 
each region from eight to four. The eigenvalue equation 
Hf = Ef entails then that the amplitudes A^_^^^{ki, fc2) 

and A^, ^, (fci , fc2) of adjacent regions i?, R' are related by 
the S'-matrices, 



(20) 



We have in general S{R, R') = S{R' , i?)"^ The S'-matrix 
corresponding to scattering between the two particles 
S'([012], [021]), taking particle 1 from the right to the 
left of particle 2, is diagonal in the pseudospin space due 
to conservation of total momentum, 



5^;:h[oi2],[02i]) = 44 '^'^^ 



sin (72^2 + iU 12 



sin (Ti fci — sin (72^2 — iU /2 
and written as a 4 x 4 matrix, it reads. 



^i2(fci,fc2) = ^([012],[021]) 



= 5(120], [210] 



/ a 

p-^ 



Vo 







/? 




(21) 



with 



sin fci — sin fc2 + ill /2 
sin fci — sinfc2 — iU /2^ 
sin fci -I- sin fc2 -I- iU /2 
sin fci + sin fc2 — iU /2 



(22a) 
(22b) 



Note that S'i2(fci, fc2)^^ = S'i2(fc2,fci) and 6*12 is unitary 
for real fci , fc2 . 

The S'-matrix describing the interaction between 
particle 1 and the potential —V reads S'io(fci) = 



12 



5([102], [012]) = S'([210], [201]). It depends on ki but 
not on k2 and has the form of 



To study the kernel of Sd further, we define 



Sio{ki 



1 + 71 71 

— 7i 1 — 7i 

1 + 71 71 

— 7i 1 — 7i 



(23) 



with 7i = —iV/{2shiki). Note that ^lo is non- unitary 
(in contrast to ^12) but preserves the current across the 
impurity: 



4(102] 



aIW2] 

— lCJ"2 



4(012] 

1(72 



4(012] 



(24) 



In an analogous way, we find for S'2o(fe2) 
^([201], [021]) = ^([120], [102]), 



^"20(^2) 



1 + 72 





72 








1+72 





72 


- 72 


1 


- 72 








- 72 





1 " 72 



(25) 



with 72 = -iV/{2smk2). Naturally, [S'io,S'2o] = 0. 

A necessary condition for the ansatz (fTO]) to work is 
its self-consistency: Starting from the amplitudes in one 
region, say [021], we can compute the amplitudes in all 
other regions using the S'-matrices. The region [120] can 
be reached from [021] on two different paths, namely 
[021] -> [012] [102] -> [120], which gives 

^[120] ^ ^2o5io^l2v4l"'^l (26) 

or via [021] [201] [210] [120], corresponding to 

A^'^^^^S^2SwS2oA^°^'^- (27) 

If an arbitrary eigenstate of Q has the Bethe form 
(dH), Eqs. ([261) and would entail the matrix iden- 
tity (Yang-Baxter equation), 

5'2o(fc2)S'io(^i)S'i2(^i,fc2) = Si2{ki,k2)Sio{ki)S2a{k2). 

Now, because [S'12, <5'iO'S'2o] 7^ 0, it is clear that the Bethe 
ansatz fails for general eigenstates. However, the matrix 

Sd = 'S'2o<5'iO'S'i2 — <5'i2'S'io<S'2o (28) 

has a two-dimensional kernel /C = {A'' ^^ A'--^^) for arbi- 
trary ki,k2,V,U and it follows 

i'^' = ^20^10^12^'^'^ = Si2SioS2oA^'\ j = 1, 2. (29) 

If the amplitude vector in region [021] is an element of 
/C, we obtain along both paths the same amplitudes in 
region [120] and the ansatz is consistent. This is also 
valid for the other possible consistency relations because 
all regions are connected by a single loop and S{R, R') = 
5(i?',i?)-i. 



= 7172(1 + a)(l - a ^). 

= 7172(1 + /3)(i-r'), 

= -(l + 72)<5//?, 
= -(1 - 72)<5/a, 
--(l + 7i)'5, 
= -(l-7i)<5/(«/3)- 



(30a) 
(30b) 
(30c) 
(30d) 
(30e) 
(30f) 



Here we have set S = {af3 — 1)71 — {a — /3)72. The 
. . . , f satisfy the relation ab + cd = ef. The matrix Sd 
reads 



Sr 



and has the characteristic polynomial. 






c 
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—a 


c 





~b 
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d 






(31) 



det(S'D - A) = (A^ 



Its kernel /C is therefore two-dimensional. With the re- 
lation ab + cd — ef we obtain two linearly independent 
basis vectors of K., 



Ad 





(32) 



The amplitudes in region [120] follow from Eq. ((29)) . 



id) 



This means A, 




(Tl ,(T2 



i(2) 



= — A[f^^__^2 3'i^d entails to- 
gether with the ansatz (|19p the relation f{xi,X2) — 
~ f{—xi, —X2), i.e. the wave function is of odd parity. 
Because the functions of the form given in Eq. p9p span 
the whole Hilbert space and are determined by a single 
region, the odd functions corresponding to elements of IC 
in region [021] span the complete odd subspace. It fol- 
lows that the projection of the Hamiltonian onto the odd 
invariant subspace can be diagonalized with the Bethe 
ansatz, the partial waves show no diffraction and the sys- 
tem is therefore integrable if confined to this subspace. 
Let us note that the exceptional case xi = X2 = has 
not to be considered separately as the wave function van- 
ishes there and the Yang-Baxter equation for two-particle 
scattering is sufficient to prove consistency of the ansatz 

m- 

The integrability of the odd subspace and the nonin- 
tegrability of the even subspace can be seen in Figs. [S] 
and [151 respectively, which depict parts of the spectral 
graph. There are apparent avoided crossings in the even 




(33) 
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Poisson 
Wigner-Dyson 




FIG. 12. (Color online) Level spacing statistics of the odd 
(blue bars) and even (red bars) subspace. The green solid line 
and the red dotted line indicate the Poisson Pp{s) = exp(— s) 
and Wigner-Dyson distribution for the orthogonal ensemble, 
Pwd{s) = exp(— 7rs^/4), respectively. The parameters are 
{V, U) = (2, 2). The lattice has 201 sites with open boundary 
conditions. Each subspace contains roughly 10000 levels. 



subspace in Fig. ^ indicating nonintegrability according 
to the level-crossing criterion 30] , while in Fig. [T5] the 
odd-parity bound state crosses both odd and even states. 
The latter is clear from parity symmetry, the former is 
due to integrability of the odd sector. This agreement 
between the numerical results on a finite lattice (which, 
however, preserves the parity symmetry) and our anal- 
ysis of the infinite lattice hints that the odd sector re- 
mains integrable on finite lattices. To look at this ques- 
tion from a different point of view, we show in Fig. [T^] 
the statistics of the level spacing in each sector. The odd 
subspace appears to be almost Poissonian, as predicted 
by the Berry- Tabor criterion for integrability [31|. The 
even subspace, on the other hand, does not show Wigner- 
Dyson statistics — characteristic for fully developed chaos 
[3^ . The statistics is close to Poissionian for large spac- 
ings as in the odd subspace, but deviates from both Pois- 
sonian and Wigner-Dyson for small spacings. Moreover, 
the numerical analysis reveals even states which have a 
close connection with associated odd states (see Sec. IVIII 
below). The intermediate behavior of the level distribu- 
tion in the nonintegrable sector may be related to a pos- 
sible "weakly diffractive" dynamics which deserves fur- 
ther investigation. For example, it would be interesting 
to check whether the classical (continuum) limit of the 
model is ergodic or not. 



C. Two Bethe form odd-parity bound states 

Having established the Bethe ansatz solvability of the 
model in the odd subspace, we now proceed to study the 
odd-parity localized states in the Bethe ansatz frame- 
work. Note that both the ground state and the molecule- 
type bound state encountered above are even and thus 
not in the Bethe form. 

To determine the conditions on ki and k2 imposed by 



the normalizability of the wave function, it is convenient 
to avoid explicit use of the symmetrizer S and go back 
to the ansatz in Eq. (fTT|) . For localized states, it is easy 
to see that among the four coefficients i?i<j<4, only one 
can be finite. The reason is that, in region I2, one of the 



two terms of 



and 



diverges to infinity while 



the other converges to zero as \xj\ — ?> 00. Therefore, their 
prefactors cannot be finite simultaneously. Without loss 
of generality, suppose Bi — 1 while B2 — — B4 ~ 
0. The (unnormalized) wave function in region I2 then 
reduces to. 



^—ik2Xi—ikiX2 



(34) 



Since in region I2, xi goes to —00 while X2 goes to +00, 
we need |e*'^^| > 1 > |e**''2| for / to vanish at infinity. 
Note that we can also define — /c2 as fci, and —ki as k2. 
To avoid counting a solution twice, we take the conven- 
tion that |e''''i+*'=2 1 < 1. From Eq. we get two inde- 
pendent equations for [/, V, ki^2, 



{U + 2i(sinA;i - sinfc2)) yl - 
{U + 2z(sin fci - sin fc2)) (1 



V 



2i sin fci 
V 



2i sin /c2 

From ((T^ and ([13]), we get furthermore, 
V . V 



0, (35a) 
: 0.(35b) 



A5 = Ae 



2i sin/ci ' 
0, = 



A, 



V 



2i sin k2 



2i sin fci 



A2 = Aa 



-1 



V 



2i sin fc2 



(36) 



In region Ii (0 < xi < 2:2), xi and a;2i = X2 — xi are 
independent variables taking values on Z+ UO. Therefore, 
we rewrite the wave function as 

fi^i X2) — Aie^^^^^^''^^^^^^^''^^'^^ -\- 

Now for the wave function to be localized, we need 
|g»fci-h»fc2| < 1 and As = 0. We thus get from dM]) 



V ^2i sin fca = eJ'"' - e 



(38) 



Eq. ((35b ) is then satisfied. For ((35k ) to be satisfied, we 
need either —V + 2i sin fci = or U + 2i (sin fci — sin fc2 ) = 
0. But the former is unacceptable, otherwise the con- 
dition |e''^i| > 1 and |e*'=i+*'=2| < i ^g^jj ]-,g satisfied 
simultaneously. Therefore, U + 2i{smki — sinfc2) — 0, or 



V ~U ^2i sinfci 



-iki 



^iki 



1 and Z2 = , we have — V 



Z2- 



(39) 



2 ' 



Defining zi = e 

—V — U = zi — z^^, and E = —Z\ — z-^" — Z2 — z^ 
with \z2\ < 1 < \zi\ < \z2\~^- Instead of studying for 
what kind of (V,U) pairs, we have solutions of (21,22) 
satisfying this condition, we do the inverse. Because V > 
by assumption, we have < Z2 < 1. Depending on the 
sign of zi, we have two cases: 
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(i) < Z2 < 1 < ^1 < z^^ We get < -V - ?7 = 
zi - < z^^ - Z2 = V. That is, -2V < U < ~V < 0. 
The energy of the wave function is 

Eb2 = - + - ^{V + Uy + i. (40) 

It is easy to prove that Eb2 < {-4, —y/V^ + 4 - 
2, —^/If^ + 16}. Consequently, this bound state is below 
all the three continuum bands and is thus not a BIC. A 
notable feature of this state is that on the line xi+X2 ~ 0, 
f{xi,X2) — and if xi + X2 > 0, f{xi,X2) < while if 
xi + X2 < 0, f{xi,X2) > 0. That is, the wave function 
has a node line xi+ X2 = 0, and is positive on one of the 
two half-planes, while negative on the other. This prop- 
erty can be readily verified from the expressions p4p and 
(EZj). 

(u) -z^^ < zi < -1 < < Z2 < 1. We get 

> -1/ - [/ = zi - 1 > -z^i + Z2 = -V. (41) 

That is, —V < U < 0. The eigenenergy of the wave 
function is 

Ebi = - + + y/{V + Uy + 4. (42) 

Now, it is easy to show that > Ebi > {U, -VV^ + 4 + 
2}, and thus Ebi falls outside of the second and third 
bands. But it can fall in the continuum band [—4, -1-4] 
to be an embedded eigenvalue. For example, for (V, U) = 
(2,-0.5), Ebi = -0.3284, which is inside the [-4,-^4] 
continuum. The condition for this state to be a BIC is 
Eb > —4. Note that this bound state exists whenever 
V < U < 0. However, only in a subset [the region de- 
noted as (ii)-BIC in Fig. [T3] of this region does it fall 
in the continuum. Its energy can be tuned continuously 
across the edge —4 by tuning the parameters across the 
dashed line in Fig. [TS] 

In Fig. [131 we have plotted the squared wave functions 
of the two bound states with a fixed V but different values 
of U. There we see that for a fixed value of V, the local- 
ized state becomes extended along the line of xi = X2 (see 
Fig.fHk) as [/ — ?> 0~, while it becomes extended along the 
lines of xi = and X2 — (see Fig. [Hh ) as U ^ —V^. 
This follows from \ziZ2 \ 1 or \zi \ 1, respectively, in 
the two limits. Similar behavior is displayed by the other 
odd-parity bound state as C/ — >■ —2V~ or —V'^, due to 
the same reason. As U crosses —V from —V" to — 
(see Figs. [T3h and fT3H). the first bound state becomes 
extended and disappears, while the second one appears 
starting from an infinite size. However, the two states 
are not continuously linked. Actually, from (15^ we see 
that —V — U plays the role of an effective defect poten- 
tial. As it changes sign, zi changes sign and the bound 
state jumps discontinuously from above the second band 
to below it. Finally, it should be stressed that in contrast 
to the power- law decay of the wave function in [l^ , here 
the wave function decays exponentially as |a;i^2| tend to 
infinity, as seen in ([Ml) and (1571) . 

Above we have shown that Ebi falls in the continuum 
band for appropriate parameters. It would be instruc- 
tive to see how this energy level moves inside the band 
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FIG. 13. (Color online) In the region labeled with (ii)-BIC, 
the state with energy Eti [see ()42|) ] is a bound state embed- 
ded in the [—4, -1-4] continuum. In the region labeled with 
(ii)-BOC, this state is a bound state below the [—4, -1-4] con- 
tinuum band, and outside of all the three continuum bands. 
The boundary between the two regions is determined by the 
condition Ebi = —4. In the region labeled with (i), the odd- 
parity bound state with energy Eb2 exists, which is below all 
the continuum bands. 
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FIG. 14. (Color online) Images of the squared wave functions 
of the two analytically solvable bound states. The value of U 
decreases monotonically from (a) to (f), with the value of V 
fixed to 2. Note that in (a)-(c), -V < 17 < while in (d)-(f), 
—2V < U < —V. These states are obtained in a finite lattice 
by exact diagonalization, with 100 sites on each side of the 
impurity and open boundary conditions. They agree with the 
analytic expressions in Eqs. ([Mil . ([551) . and ([571l . 



as the parameters vary. Since the BIC is an eigenstate 
of the integrable sector, it is expected that Ebi will cross 
the levels of the extended states without hybridization. 
This is indeed the case. In Fig. [151 we have shown a 
portion of the spectral graph of the model. We see that 
the BIC level (the diagonal one) is the only level that 
changes significantly in the interval of U and it crosses 
all other (even or odd) levels leading to degeneracies at 
the crossing points. For even extended states it is a con- 
sequence of the parity-symmetry, while for odd states it 
is expected from integrability of the odd sector [soj . 
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FIG. 15. (Color online) A small portion of the spectrum of 
the model with V = 2. The lattice has 201 sites with open 
boundary conditions. In total, 10 levels are shown. The di- 
agonal one corresponds to the bound state in the continuum 
(BIG). The horizontal lines correspond to extended states in 
the [—4, +4] continuum band. Some of them are even, some 
are odd. The BIG level, which is odd, crosses both the even 
and odd ones. 
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(a) V = 1.5, U = -2 




FIG. 16. (Color online) The two analytically solvable, 
odd-parity bound states and their accompanying even-parity 
bound states (each pair is indicated by an arrow) revealed by 
using the lattice expansion algorithm. The values of (V^, U) 
are shown in each panel. The upper panel corresponds to the 
case of —2V < U < —V, while the lower panel corresponds 
to the case of ~V < U < 0. The smaller and larger lattices 
have 51 and 71 sites, respectively. 



VII. TWO ACCOMPANYING BOUND STATES 

In the proceeding section, we have solved two odd- 
parity bound states analytically, each of which appears 
in an appropriate region. They appear in Figs. and 
Numerically, it is observed that often each of the 
two states is accompanied by an even-parity bound state, 
which is close to it in energy. Actually, in Fig. [Sjj, one of 
the excited bound states has odd-parity, while the other 
has even-parity. In Fig. we see only the odd-parity 
one (the BIC). However, if its energy is tuned outside 
the [—4, +4] continuum band, it can have an accompa- 
nying even-parity bound state. In Fig. I16[ we show two 
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FIG. 17. (Color online) (a) and (b): Images of the squared 
wave functions of the two analytically solvable bound states 
in Fig. 1161 and (d) and (e): their accompanying even-parity 
bound states. The values of {V, U) are the same as in Fig. 1161 
i.e., {V,U) = (1.5,-2) in (a) and (d), or {V,U) = (5.7,-5.3) 
in (b) and (e). These states are obtained in a lattice with 201 
sites and open boundary conditions by exact diagonalization. 
In (c) and (f), the squared wave function of the sum or the 
difference of the two bound states in (b) and (e) are shown, 
respectively. 



cases when the two odd-parity bound states each has an 
accompanying even-parity bound state. 

Unlike the odd-parity bound states, the even-parity 
bound states have not the simple Bethe form. This makes 
them much more difficult to study and understand. How- 
ever, thorough numerical exploration reveals that they 
are closely related to their odd-parity companions. First, 
in each pair, the two states are ordered in a definite way 
and the gap between them is always much smaller than 
unity. Specifically, the odd-parity state with energy Eb2 
(if it exists) is always the first excited state, while its 
companion is (if it exists) the second excited one. For 
the odd-parity state with energy Ehi^ its companion (if 
it exists) is always lower in energy to it. Second, as il- 
lustrated in Fig. [ITl the wave functions of the two bound 
states in each pair are similar to each other in profile. 
More precisely, if we denote the wave functions of the 
even and odd states in each pair as fe and /o, respec- 
tively, then fe + fo and fe — fo are well localized on one 
of the two half-planes xi -t- 2:2 ^ 0. This is demonstrated 
in Figs. [T7h and [TTT . These two facts remind us of the 
double-well potential, where we have pairs of even and 
odd eigenstates, which are nearly degenerate in energy 
and can be expressed as linear combinations of left- and 
right- localized modes. It indicates that the wave func- 
tion of the even-parity bound state, while not exactly in 
the Bethe form, may be nevertheless close to it, a numer- 
ical finding which should be further studied. 

We have tried to map out the regions where the two ac- 
companying bound states appear. The results are shown 
in Fig. [181 In Fig- ITSk . we see that if U is too close to 
the edges of the region —2V < U < —V, the even-parity 



16 




V V 

FIG. 18. (Color online) (a) In the region enclosed by the 
red lines, an accompanying even-parity bound state to the 
analytically solvable odd-parity bound state with energy Eb2 
exists. The green dotted lines enclose the region where the 
simple variational approach predicts the existence of such an 
accompanying bound state. The blue dashed lines enclose 
the region —2V < U < —V where the Eb2 bound state exists, 
(b) In the region enclosed by the red lines, an accompany- 
ing even-parity bound state to the analytically solvable odd- 
parity bound state with energy Eti exists. The Eti bound 
state exists below the straight blue dashed line and moves 
into the continuum below the curved line (comp. Fig. I13p . 

bound state disappears. The reason is that in these two 
limits, i?b2 gets close to the lower edges of the second or 
third band [see Eq. (|40p ]. The even-parity one would get 
even closer or just inside the bands and thus becomes 
unstable. As for the second even-parity bound state, its 
appearance region is approximately the region where the 
Eti odd-parity bound state is a BOC. As Ebi gets too 
close to the upper edge of the second band {U -> —V^), 
or too deep inside the [—4, +4] band, the even-parity 
companion becomes unstable and disappears. 

Here, we would like to remark that for the odd-parity 
bound state with energy Ei,2 , the presence of an accompa- 
nying even-parity bound state is expected in some region. 
As mentioned in Sec. IVIB . the wave function of the odd- 
parity bound state fb2 vanishes on the line xi -f a;2 = 
and is completely positive on one of the planes xi+X2 ^ 
while completely negative on the other. Now construct 
an even-parity state 

r _ j +.fb2{xi,X2), a;i+.T2>0, 

[~Jb2[Xl,X2), Xi+X2<0. 

By Eq. ®, we have then {ftr\H\ftr) = {fb2\H\fb2) = 
Eii2, which is lower than the lower edges of all the 
three bands. Then construct another even-parity state 
ftr ~ (1 ~ \g}{9\)ftr- This state is even and orthogonal 
to both the ground state \g) and the odd-parity state /f,2. 
If its energy is lower than the lower edges of the three con- 
tinuum bands, then there must exist another even-parity 
bound state. This is likely since if -2V <U < -V < 0, 
the ground state \g) localizes sharply at the origin where 
ftr is zero and thus ftr is actually very close to ftr- Ac- 
cording to this criterion, we can prove that in the re- 



gion enclosed by the green dotted lines in Fig. [T5k . there 
should be an even-parity bound state. 



VIII. CONCLUSIONS AND OUTLOOK 

We have investigated the bound states in a system of 
two interacting bosons in a one-dimensional lattice with 
a defect. The aim was to study a problem known from 
atomic physics on a lattice. As the simplest one among 
its type, our model is both analytically and numerically 
tractable. Somewhat unexpectedly, this simple model 
reveals rich physics. 

First, we find that the repulsion between the two 
bosons plays a dual role. On one hand, as intuitively 
expected, a sufficiently large U can delocalize the ground 
state. On the other hand, a sufficiently large U can cre- 
ate a bound state and stabilize it [331 ■ Quantitatively, 
there exist two critical values C/ci,2 (0 < f/ci < Uc2) of 
U. For U > Uci, the ground state is no longer local- 
ized. However, for U > Uc2, a molecule- type bound state 
appears and becomes more and more localized as U in- 
creases. The large U can in fact bind the two bosons 
into a molecule, which hops on the lattice with an effec- 
tive amplitude cx \/U, and can thus be trapped by the 
defect. This localization becomes sharper with increasing 
U because the effective hopping decreases as 1/U . 

An interesting point is that the two transitions have 
distinct natures. For the former, the ground state energy 
and wave function merge into the continuum continu- 
ously as U increases towards Uci- Near Uci, the wave 
function can be well approximated by a modified Chan- 
drasekhar wave function and the value of Uci can be ac- 
curately predicted in this way. In contrast, for the latter, 
the bound state energy merges into the continuum by 
a sharp anti-crossing as U decreases towards Uc2, which 
means that the bound state becomes extended almost 
abruptly. This bound state is a hound state at threshold 
, 13-17] ■ which has a finite size even at the threshold of 
the continuum. 

Second, we have gained a better understanding of the 
analytical structure of the model. Though the model was 
believed to be nonintegrable [2C|, we have seen that the 
odd-parity subspace is integrable. This makes it possible 
to write down the wave functions, in particular, the two 
bound states, explicitly in the Bethe form. Remarkably, 
one of the two turns out to be a hound state emhedded 
in the continuum (BIG). It is the simplest BIG we know 
so far, in terms of the energy and wave function; and 
most importantly, it arises in a non-engineered model. 
Therefore, in this simple model, there can be two types 
of exotic bound states: one remains finite at the threshold 
of the continuum; the other resides inside the continuum. 

Eigenstates in the even-parity subspace do not have 
the Bethe form. However, the fact that the level spac- 
ing statistics in this subspace interpolates between the 
Poisson distribution and the Wigner-Dyson distribution, 
hints at a weakly diffractive dynamics even in this non- 
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integrable sector, which should be investigated further in 
the future. 

Third, besides the four types of bound states men- 
tioned above, we have found numericahy that in some re- 
gions of the parameter space the odd-parity bound states 
may possess accompanying, even-parity bound states, i.e. 
there are two bound states of different parity which are 
close in energy and have similar wave functions. The situ- 
ation resembles the nearly degenerate pairs of eigenstates 
in a double-well potential. Again, this phenomenon im- 
plies that there exists some relation yet to be discov- 
ered between the diffractive even sector and the non- 
diffractive odd sector. 

In this paper, we have focused on the properties of 
the bound states, which is mainly a static problem. In 
the future, it is worthwhile to study the dynamics of the 
model [U, |35|. In fact, the analogy of the model with 
some atomic physics problems can be generalized to the 
dynamics. For example, consider the following: Initially 
one boson is trapped in the defect mode and the other 
boson comes in from faraway. What would the system 
evolve into? There are many possibilities. First, the 
injecting boson can be reflected or transmitted without 
exciting the initially trapped boson. Second, the initially 
trapped boson may get excited and both bosons become 
free and move independently. Third, the two bosons can 
go away as a molecule. It is a nontrivial task to find the 
probability for each of these possibilities. The problem 
obviously bears analogy to the well-known problem of 
electron scattering off hydrogen atoms . 
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Appendix A: Variational proof of Uci > V 

We argue that if [/ < V, the ground state is a bound 
state. To this end, we just need to prove that the ground 
state energy is below all the three continuum bands, since 
any eigenvalue outside of the continuum bands should be 
a bound state according to our argument in Sec. |lTl 

Consider a direct product variational wave function 
f{xi,X2) = 0i(xi)(/()2(x2), with 0i_2 normalized to unity. 
Its energy E = {f\H\f) is 

E - -I- {h\h\h) + uY,4>\{x)cl>l{x). (Al) 

X 

Here h is the single particle Hamiltonian, i.e., the sum 
of the hopping and the potential part of the full Hamil- 
tonian. Now it is clear that particle 1 sees an effective 
potential Wi{x) = —VSxfi + U4>\{x), which depends on 



the wave function of particle 2, while particle 2 sees an 
effective potential W2{x) = —VS^.o + U(j>i{x), which in 
turn depends on the wave function of particle 1. Re- 
gardless of the concrete forms of (/)i.2, we have always 
X]m^i("^) — Z]m^2(w) = U —V. Variationally we 
know that if U < V, then both Wi and W2 admit local- 
ized ground states. Now choose (pi as the ground state of 
h with eigenvalue —y/V^ + 4, and then choose 02 as the 
ground state, a localized state, in the effective potential 
W2- We have then. 



E = {q^i\h\(j)i} + {cP2\h + W2\(j)2) < -VV^+l~2. (A2) 

We can decompose the direct product wave function as 
the superposition of a symmetric (bosonic) and an anti- 
symmetric (fcrmionic) function, f{xi,X2) = fs{xi,X2) + 
faixi,X2), with 

fs{xi,X2) = ^ (0i(a;i)02(a;2) + <j)2ixi)(j)i{x2)) , (A3) 

fa{xi,X2) = ^ (01 (Xi) 02 (2:2) - 4'2{xi)(f>l{x2)) . (A4) 

Note that fs and fa are not normalized. Actually, 1 = 
(/I/) = ifslfs) + (fa\fa)- The even and odd functions do 
not mix under the action of H . Therefore, 

(/l^l/> = (/.|/s)%^ + (/a|/a)%^. (A5) 



(/sl/s 



(/a I /a) 



We note that for a fcrmionic wave function as /q, the 
interaction is non-effective. Therefore, by Pauli principle, 
we have 



{fa\H\fa) 



{fa\fa 

From (|X2|) and (KE\i . we have 



> -y/V^+4:-2. 



{fs\H\fs 
ifslfs) 



(A6) 



(A7) 



Therefore, we have constructed a localized bosonic state 
fs, which is below all the three bands. Thus, the ground 
state is a localized state iiU < V. We have thus obtained 
a lower bound for Ud, i.e., Ud > V. 



Appendix B: Lower limit of Uci by Chandrasekhar 
trial wave functions 

We employ the Chandrasekhar trial wave function 

fixi,X2) = e^"il^il-"2l^2l -i-e-^^l^'il^^il^^l, (Bl) 

with two variational parameters satisfying < ai < a2- 
Note that the wave function is in the form of fs [see 
Eq. (jA3p ] in App. |21 but with the orbits 0i^2 specified. 
We have for the norm of the wave function 



(/I/) 



tanh ai tanh 02 tanh a tanh a ' 



(B2) 
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with a = (ai + Q!2)/2 and 5 = {a2 — We have also 

4 4 

{f\H\f) = 



sinh ai tanh a2 sinh a2 tanh ai 
8 cosh (5 4C/ 



sinh a tanh a tanh(ai + 012) 



'2V 



1 



1 



tanh ai tanh a2 tanh a 



.(B3) 



The energy of the trial wave function is then E 
{f\H\f) /{f\f)^ (a + b)/ic + d), with 



4 



2V 



sinh ai tanh a2 sinh 02 tanh ai tanh ai 
8 cosh J 



sinh a tanh a tanh(Q;i +0:2) 



\ tanh a2 tanh a ^ 



tanhai tanh a2 



tanh a tanh < 



,(B4) 

(B5) 
(B6) 



Since c,d > 0, we have E > min (a/c, fc/d). However, 
it is easy to prove that a/c > Ef,dge2 by noting that 
sinh(ai) > tanh(ai). Therefore, we need h/d < Eedge2 
to make E < Eedge2 possible. Numerical evidence indi- 
cates that in this variational framework, at the critical 
point oi U = Uci, ai = and ^2 = Hhi^ + + 4)]. 
That is, one boson is delocalized while the other resides 
in the single-particle localized mode corresponding to the 
defect potential. 

The critical value Uci is then determined by the condi- 
tion b/d = Eedge2 with ai^2 taking the values above. It is 
straightforward to obtain the limiting behaviors of Uci '■ 




(B7) 



It is an improvement over the variational approximation 
in App. in both limits. However, it is still inferior to 
the modified Chandrasekhar variational approximation. 
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